Could the Casimir Effect explain the Energetics of High- Temperature 

Superconductors? 
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It is proposed that the Casimir energy of the electromagnetic field in between the layers of cuprate 
superconductors could significantly contribute to the energy balance in these materials. Since the 
Casimir energy is strongly dependent on the distance between the layers a corresponding dependence 
is predicted for the superconducting transition temperatures. Comparison with the experimental 
data on the transition temperatures of epitaxial superlattices of alternating layers of YBCO and 
PrBCO of varying thickness shows that these are well reproduced. 
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PACS numbers: 74.72.-h, 74.20.-z,12.20.-m 

The question has remained open how Cooper pairs 
can be stable at around lOOK where some high tempera- 
ture superconductors (HTSCs) are still superconducting. 
In particular, the phonon-mediated attractive electron- 
electron interaction of BCS theory is known to be too 
weak at these temperatures, see, e.g., 

For a new approach to this problem, let us reconsider 
a feature that HTSCs have in common, namely paral- 
lel superconducting layers which are separated by layers 
of essentially insulating material. Since in between any 
two conducting surfaces there occurs a Casimir effect, see 
[H, the effect also occurs between the parallel supercon- 
ducting layers in HTSCs, as was first pointed out in 
Before estimating its significance in HTSCs, let us briefly 
review the Casimir effect's underlying mechanism. 

We begin by recalling that the ground state energy of 
a quantum harmonic oscillator is not necessarily fixed. 
It can be lowered, for example, by coupling the oscilla- 
tor to degrees of freedom which decrease the amplitude 
of its zero-point fluctuations. This happens, e.g., when 
these extra degrees of freedom effectively increase the 
oscillator's mass. An example would be the vibrational 
oscillator of a diatomic molecule that captures a neutron. 
Similarly, the ground state energy of the electromagnetic 
field can be lowered by suppressing the zero point fluc- 
tuations of some of its electromagnetic field modes. As 
Casimir showed, this happens, for example, when the 
electromagnetic field couples to conducting charge car- 
riers which are confined to two parallel plates, in which 
case certain modes in the c-direction become suppressed. 
In the case of HTSCs, as the temperature is lowered be- 
low Tc, superconducting charge carriers form in paral- 
lel layers and they therefore lower the electromagnetic 
ground state energy. Our aim here is to estimate if this 
lowering of the energy might be sufficient to make the 
formation of Cooper pairs energetically favorable at tem- 
peratures as high as lOOK. 

The derivation of the Casimir effect for two conduct- 
ing objects from first principles requires the calculation 



of the ground state energy of the quantum system that 
consists of both, the electromagnetic field and the charge 
carriers. This calculation is hard and it is instructive 
to consider first the simplified case where the two ob- 
jects are what may be called ideal conductors, namely 
conductors whose conductivity and MeiBner effect ex- 
pel electromagnetic fields of all wavelengths with van- 
ishing penetration depths. In this case, the charge de- 
grees of freedom are effectively integrated out and one 
is left with the noninteracting quantum field theory of 
electromagnetism together with the boundary condition 
that the electromagnetic fields vanish at the objects' sur- 
faces. Only those electromagnetic modes which obey this 
boundary condition contribute their ground state energy. 
In this way, the total ground state energy of the system 
depends on the objects' (shape and) distance, a. The 
distance dependence of the energy implies a force, the 
Casimir force. These calculations require renormaliza- 
tion but are straightforward for simple geometries such 
as two parallel plates. The Casimir force between two 
ideally conducting, neutral and parallel plates of large 
area A and distance a is given by: 
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Corrections that take into account the finite conductiv- 
ity of real metals have been calculated for geometries 
such as parallel plates and a plate and a sphere, along 
with corrections for finite surface roughness and finite 
temperature, see |^. Recent experiments measured the 
Casimir force between a metallic plate and sphere down 
to distances around lOOnm, confirming the theoretical 
predictions to a precision of 0.5% 6J. 

From Eq^ the ground state energy of the quantum 
system consisting of the electromagnetic field and two 
ideally conducting, neutral and parallel plates reads: 
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The integration constant is chosen such that the energy 
of the system vanishes at infinite plate separation, or 
equivalently, in the absence of conducting plates. We see 
that a system of two conducting parallel plates at a dis- 
tance a is energetically lower than the same system of 
two parallel plates at distance a if the two plates are in- 
sulators. Indeed, if two parallel insulating plates possess 
a microscopic mechanism which allows them to create 
ideally conducting charge carriers at an energy expense 
of no more than E(a) (and correspondingly less at fi- 
nite temperature), then the two plates are energetically 
driven to use this microscopic mechanism to create ide- 
ally conducting charge carriers. 

This scenario may apply to HTSCs. They possess par- 
allel layers which are initially insulating but can become 
superconducting at low enough temperature, i.e., there 
exists some microscopic mechanism that allows these lay- 
ers to create superconducting charge carriers. The for- 
mation of the Cooper pairs leads not to ideal conduc- 
tivity but to superconductivity, which in turn leads to a 
partial suppression of the fluctuations of certain electro- 
magnetic modes. This leads, therefore, to some lowering 
of the vacuum energy of the electromagnetic field. If this 
lowering of the electromagnetic zero-point energy at the 
onset of superconductivity is large enough then the Cu-0 
layers' initially non-superconducting charge carriers are 
indeed energetically driven to utilize the available micro- 
scopic mechanism in order to become Cooper pairs. We 
notice here the importance of the non-superconducting 
state being an insulating state. This is because the size 
of the change of the electromagnetic ground state energy 
is crucial and it depends on the difference in conductivity 
in the two phases. 

In this scenario, it would be the very effects of su- 
perconductivity which enable and stabilize superconduc- 
tivity. Cooper pairs would derive their stability collec- 
tively, across layers. Namely, Cooper pairs would be sta- 
ble because if sufficiently many of the Cooper pairs on 
opposing layers were to break up then the suppression 
of some electromagnetic modes' zero-point fluctuations 
would cease and their electromagnetic ground state en- 
ergy would have to go back up to its unsuppressed level. 

While a Casimir effect must occur in HTSCs it is, a 
priori, not clear if it is indeed strong enough to lead to 
this scenario, especially because the superconducting Cu- 
O layers are less efficient at suppressing electromagnetic 
modes than ideal conductors would be. Let us, there- 
fore, consider the well-studied material YBa2Cu307^x 
(YBCO) which becomes superconducting at around 92K. 
The crystallographic unit cell contains two copper oxide 
layers at a distance of ab ~ 0.39nm and neighboring such 
bi-layers are separated by a layer of essentially noncon- 
ducting material of width ~ l.llnm. The area density 
of superconducting charge carriers on each individual Cu- 
O layer may reach on the order of rig « 10^''/cto^. 

For our purposes, the case of YBCO is of particu- 



lar interest because of the availability of experimental 
data, 0, 0, on epitaxial superlattices in which slabs 
of YBCO alternate with slabs of insulating material, 
namely PrBa2Cu30T-x(Pi^CO). For example, in the 
experiments reported in |g] , the authors varied the thick- 
ness of the YBCO slabs from 7\/ = 1 to = 8 unit 
cells and the thickness, Um, of the PrBCO slabs from 
iV = 1 to = 16 unit cells, i.e., in the range Um = 2nm 
to 20nm. The superconducting transition temperature 
was measured as a function, Tc{N,am), of N and a,„. 
The results are summarized in Fig. 3 of 0. Since Tc is 
known to be a function of the layer separation and since 
the Casimir effect is sensitive to varying layer separations 
those data provide a good testing ground for the present 
ansatz. 

Let us initially consider the simplified case in which 
the Cu-0 layers in their superconducting state are taken 
to be ideal conductors in the above sense, which means 
that they cause electromagnetic fields of all wavelengths 
to vanish on their surface. We also assume that the Cu-0 
layers are of negligible thickness and separated by vac- 
uum. In this case, the charge degrees of freedom are 
effectively integrated out and we are left with the free 
quantum field theory of electromagnetic fields with the 
boundary condition that these fields vanish on each of 
the superconducting Cu-0 layers. The ground state en- 
ergy of the electromagnetic field between any two ideally 
conducting layers of distance a is then lowered by the 
amount given in Eq|21 and we can apply this result to all 
the inter-layer distances that occur in the superlattice. 

Concretely, each period of the superlattice contains a 
slab of M unit cells of YBCO (each cell with one bi-layer 
of Cu-0), followed by a slab of insulating PrBCO of thick- 
ness am- Therefore, in each period of the superlattice, the 
case that two neighboring superconducting Cu-0 layers 
are separated by the distance ab occurs M times. The 
case that two neighboring superconducting Cu-0 layers 
are separated by the distance occurs (M — 1) times 
and finally the case that two neighboring superconduct- 
ing Cu-0 layers are separated by the distance -I- 
occurs once per period of the superlattice. The reduction 
of the ground state energy of the combined electromag- 
netic and charge carrier system, within the volume given 
by one period of the superlattice times an area A in the 
ab plane, is therefore given by: 
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Since one period of the superlattice contains 2M layers 
of Cu-0, this energy reduction is shared by 2MnsA su- 
perconducting charge carriers, yielding the gap energy: 
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In HTSCs, the value of the variable rj, 
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is thought to be around or somewhat larger than the BCS 
value of T] K, 3.5. We obtain for the temperature Tc of 
the superconducting transition: 
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As Fig^ shows, the so-predicted transition temperature 
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FIG. 1: Tc for the ideal Casimir effect and rj — 3.5, — 
5x lO^Vcm^ The curves are, from bottom to top, for YBCO 
layer thicknesses M = 1, 2, 3, 4, 8, oo (the latter is the case of 
pure YBCO). Notice the temperature scale of kilo Kelvin. 

curves are qualitatively of the right shape when compared 
with Fig. 3 of 0. However, the transition temperatures 
are four orders of magnitude too high! The onset of 
ideal conductivity on these Cu-0 planes would provide 
far more Casimir energy than is needed to explain the 
Cooper pair binding energy. But the predicted Casimir 
effect for real Cu-0 layers is of course weaker. 

The expression Eq|21 for the Casimir energy function 
must be corrected, for example, to take into account that 
the space between the Cu-0 layers is dielectric rather 
than a vacuum. It is also necessary to correct for the 
non-flatness of the Cu-0 layers and for finite temperature 
effects. Most importantly, however, real superconducting 
Cu-0 layers are not ideal conductors in the sense above. 
Real Cu-0 layers merely suppress the penetration of elec- 
tromagnetic fields and they do so less and less efficiently 
the shorter the wavelength. In particular, modes whose 
wavelengths in the c direction are shorter than a certain 
length scale contribute to the ground state energy essen- 
tially the same amount whether the Cu-0 planes are in 
their insulating or in their superconducting state. These 
modes do not contribute to the Casimir effect, i.e., to 
the lowering of the ground state energy at the onset of 
superconductivity in Cu-0 layers. 



^ for a ^ 0, 
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Thus, while the Casimir energy, E{a), of the idealized 
case given in Eq|2| diverges as E{a) cx a 
the corrected Casimir energy function Eco, 
ten at some critical length, say Uc- For a ^ 0, it is to 
be expected that £'corr(a) because the volume in 
which the lowering of the density of the electromagnetic 
ground state energy occurs, namely the volume between 
the plates, goes to zero as a ^ 0. 

The critical length, Oc, should be roughly in the range 
10°...10^nm. Namely, while Oc cannot be smaller than 
the coherence length in the c direction, ac also should not 
be much larger than about lOOnm. The latter bound is 
suggested by the fact that measurements of the Casimir 
effect between non-superconducting metal objects have 
shown that the Casimir force persists even there, without 
changing sign, at least down to values of a « lOOnm, see 
[^. Let us write the corrected expression _Bcorr(a) as: 
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Here, / is a cutoff function of which we know that it 
obeys /(a) probably at least as fast as for a <C Cc, 
and we make the simplifying assumption that /(a) « z 
for a ^ ac with z — 1. It is important to calculate / 
for concrete materials but this is hard as it involves the 
renormalization of relativistic quantum electrodynamics 
in a highly nontrivial background. For now, our aim is to 
establish only roughly the potential role of the Casimir 
effect in HTSCs and for this limited purpose the pre- 
cise behavior of / should not matter. Also, while the 
value of z is likely smaller than 1, any value of z is 
straightforward to accommodate in our subsequent cal- 
culations. Below, we will therefore consider the simple 
example /i(a) = e^°'^^°-, and for comparison also the ex- 
ample /2(a) = e~^f°^^ which describes a slightly softer 
cutoff. In a further simplification, we treat the Casimir 
effects of neighboring pairs of layers as independent so 
that we can then use / to calculate £'iorr"°'^'' as before. 
We obtain for the critical temperature: 



^ A if (at) + (1 - l/M)/(aO + /(a„)/m) 

= 1440 ks 
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At this point, both ry and rig are still free parameters, but 
we notice that only their product enters in Eq|Sl The 
question is, therefore, if, given a cutoff function, all six 
experimental curves reported in Fig. 3 of ^] can be repro- 
duced at once when appropriately choosing the two free 
parameters ac and rj ■ Us- That these six curves can in 
fact be reproduced surprisingly well is shown in Figs l2l3l 
The comparatively large magnitude of the ideal Casimir 
effect (i.e. the case ac — 0), see Fig^ shows that, in 
principle, there is in this scenario no energetic barrier to 
raising Tc to significantly higher temperatures. The chal- 
lenge is to lower Oc, i.e., to find layered materials, such as 
cuprates, for which the Casimir effect is as large as pos- 



4 




oH — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I 

5 10 15 20 

PrBCO layer thickness [nm] 

FIG. 2: Tc predicted for corrected Casimir effect with the 
cutoff function fi(a) = e'^""^^ and setting ac = 2.7nm, t] = 
9.5, = 9 X The curves are, from bottom to top, 

for YBCO slab thicknesses of M = 1, 2, 3, 4, 8, oo. 
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FIG. 3: Tc predicted for corrected Casimir effect with the 

cutoff function /2(a) = e~^/~°^^ and setting ac = 45nm, 77 = 
0.0, Us = 5 X lO^Vcm.^ The curves are, from bottom to top, 
for YBCO slab thicknesses of M = 1, 2, 3, 4, 8, oo. 



sible. To this end, the material must normally be a very 
good insulator while in the superconducting state the su- 
perconducting layers must be able to efficiently suppress 
electromagnetic fields down to an as small as possible 
wavelength, a^- 

Of course, Tc could also be increased by reducing the 
layer separations and Oj. For example, Tc should gen- 
erally increase with pressure in the c-direction, except for 
situations where and ai become significantly smaller 
than Oc, in which case the absolute value of the Casimir 
energy drops and therefore should drop as well. This 
might occur, e.g., when the number of Cu-0 layers per 



unit cell is increased. A detailed comparison with data, 
see, e.g., 0, is in progress. 

Depending on / and the relative sizes of a^, ai and 
ac, the crystal should also exhibit a finite expansion or 
contraction in the c direction as it is cooled from just 
above to just below Tc, due to the onset of the Casimir 
force between Cu-0 layers. E.g., YBCO's elastic mod- 
ulus oiE lO^GPa yields positive AL/L « 0(10-'') 
to 0(10^^) for the case of / and Oc as in Figs. 2, 3. This 
effect could be measurable. 

In addition, seemingly paradoxical situations are to 
be expected, where a lower charge carrier density leads 
to a higher Tc, or to phenomena related to a pseudo- 
gap. This is because when there are fewer charge car- 
riers, the Casimir energy can be shared among fewer of 
them, thereby increasing the Cooper pair binding energy 
available to each pair. This is counteracted by the fact 
that the Casimir energy only materializes as long as the 
superconducting charge carrier density maintains a sup- 
pression of the electromagnetic field in the Cu-0 planes 
that is sufficient to yield a Casimir effect. 

We close with a remark on superconducting carbon 
nanotubes (CN), see 0. From the above, it is plausible 
that, due to the Casimir effect, Tc may be higher for 
multi-walled than for single-walled CNs. Indeed, very 
recently, a 30-fold increase in Tc for multi-walled CNs 
has been reported [l^ . 

In conclusion, we have shown that the Casimir effect 
could play a significant role in HTSCs, possibly even ac- 
counting for the overall energetic stability of the Cooper 
pairs. The question as to the detailed microscopic bind- 
ing mechanism for Cooper pairs would still remain largely 
open, however. We could only say that a) the stability 
of Cooper pairs would be a collective phenomenon as it 
involves mode suppression across layers and b) that since 
the Casimir effect concerns the energy stored in electro- 
magnetic modes, the relativistic retardation of the elec- 
tromagnetic interaction would play a role in HTSCs. Our 
calculations are simple enough to be qualitatively robust. 
Quantitatively, however, the validity of our results hinges 
on the assumptions made about the cutoff function / and 
on the assumption that the Casimir effects of neighbor- 
ing pairs of layers can be treated as independent, i.e., 
that these Casimir energies are additive. These points 
will need to be investigated from first principles. 

Acknowledgement: This work has been supported by 
CFI, OIT, PREA and the Canada Research Chairs Pro- 
gram of the National Science and Engineering Research 
Council of Canada. 



[1] M. Tinkham, Introduction to Superconductivity, Dover 
Publications, Mineola, N.Y. (1996), P.W. Anderson, 
The theory of superconductivity in the High-Tc Cuprates, 



5 



Princeton University Press, Princeton, N.J. (1997), Y. 

Yanase et al, Phys. Rept., 387, 1 (2003) 
[2] J.R. Schriefler, Theory of Superconductivity, Benjamin 

Inc., New York, N.Y. (1964) 
[3] H.B.G. Casimir, D. Polder, Phys. Rev. 73, 360 (1948) 
[4] A. Kempf, talk given at the 10th Marcel Gross- 

mann Meeting, Rio de Janeiro, Brazil, 20-26.7.2003, 

gr-qc/0403112 

[5] G. Barton, N.S.J. Fawsett, Phys. Rept. 170, 1 (1988), 
V.B. Bezerra, G.L. Klimchitskaya, V.M. Mostepanenko, 
Phys. Rev. A62, 014102 (2000), M. Bordag, U. Mo- 
hideen, V.M. Mostepanenko, Phys. Rept. 353, 1 (2001), 
Y Aghababaie, CP. Burgess, Phys. Rev. D70, 085003 
(2004), S.K. Lamoreaux, Rep. Prog. Phys., 68, 201 
(2005) 



[6] R.S. Decca et al.. Annals Phys. 318 37 (2005), 

quant-ph/0503105 U. Mohideen, A. Roy, Phys. Rev. 

Lett., 81, 4549 (1998) 
[7] D.H. Lowndes, D.P. Norton, J.D. Budai, Phys. Rev. Lett. 

65, 1160 (1990) 
[8] Q. Li, X.X. Xi, X.D. Wu, A. Inarn, S. Vadlamannati, 

W.L. McLean, Phys. Rev. Lett. 64, 3086 (1990) 
[9] E.V.L. de Mello et al., Phys. Rev. B 66, 092504 (2002), 

C. Acha et al., Phys. Rev. B 57, R5630R5633 (1998) 
[10] S. Chakravarty, H.-Y. Kee, K. Volker, Nature, 428, 53 

(2004) 

[11] Z.K. Tang et. al.. Science, 292 2462 (2001), M. Kodak 

et. al., Phys. Rev. Lett. 86 2416 (2001) 
[12] I. Takesue, et. al., Phys. Rev. Lett. 96, 057001 (2006) 



